Generalized Clausius inequality for nonequilibrium quantum processes 
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We show that the nonequihbrium entropy production for a driven quantum system is larger than 
the Bures length, the geometric distance between its actual state and the corresponding equilibrium 
state. This universal lower bound generaUzes the Clausius inequaUty to arbitrary nonequilibrium 
processes beyond linear response. We further derive a fundamental upper bound for the quantum 
entropy production rate and discuss its connection to the Bremermann-Bekenstein bound. 
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All real niacroseopie processes are irreversible. In 
thermodynamics, irreversibility is quantified by means 
of the entropy S: For any state transformation, the 
variation of entropy is written as /S.S = AS'i-g + AS'ii-, 
where AS'rc = Q/T is the entropy change associated 
with reversible (equilibrium) processes [T]. The Clau- 
sius inequality, AS'ir > 0, on the other hand, provides 
a fundamental characterization of irreversible (nonequi- 
librium) phenomena by specifying a lower bound for the 
irreversible entropy change; this lower bound (zero) is 
trivially independent of how far from equilibrium a pro- 
cess operates. In many cases of interest, however, hav- 
ing a sharper, transformation dependent lower bound is 
essential. A case in point is the optimization of the per- 
formance of real thermodynamic processes that occur in 
finite time [31 13] . For classical, near-equilibrium trans- 
formations, such a lower bound has been derived using 
a geometric approach to thermodynamics |4]: The in- 
finitesimal irreversible entropy production is given by the 
Riemannian distance between initial (equilibrium) and 
final (nonequilibrium) states, dS^^ > M^ /2 [315]; this ex- 
pression is obtained by a second-order expansion around 
equilibrium and is therefore restricted to the linear re- 
sponse regime. The thermodynamic length ^ measures 
the number of distinguishable states between initial and 
final probability distributions, po and p-n and is identi- 
cal to Wootters' statistical distance between wave vectors 
(pure states) in Hilbcrt space j^ ; it is explicitly given by 



HpoiPt) 



arccos 



dx y^po (x) pr (x) 



(1) 



and hence measures the angle in state space separating 
the two probability distributions po and Pt- 

In this paper, we extend the above results to quan- 
tum, far from equilibrium transformations. Specifically, 
we generalize the familiar Clausius inequality by deriv- 
ing a universal lower bound for the irreversible entropy 
production valid for arbitrary quantum processes, using 
tools from quantum information theory. We moreover 
show that, in contrast to classical nonequilibrium physics, 
there exists a maximum quantum entropy production 
rate. Our work is motivated by recent experiments on 
driven cold-atom gases that for the first time allow the di- 



rect investigation of the nonequilibrium dynamics of iso- 
lated many-particle quantum systems beyond the near- 
equilibrium linear response regime [101 111] , A precise 
characterization of the quantum nonequilibrium entropy 
production in this unchartered domain appears therefore 
necessary. In the following, we start from a recently de- 
rived formula for the total work done on a closed quan- 
tum system [12] to obtain a microscopic expression for 
the irreversible entropy production, valid far from equi- 
librium. We show that the latter is bounded from below 
by the Bures length [T31[I1], a quantum generalization for 
mixed states of Wootters' statistical distance. The Bures 
length is closely related to the quantum fidelity, a central 
measure of quantum information theory |15j . We further 
demonstrate that the quantum entropy production rate 
a is bounded from above by a quantity that also depends 
on the Bures length. This fundamental limit on the en- 
tropy variation rate is of purely quantum origin and is 
connected to the energy-time uncertainty relation. Re- 
markably, we show that the maximum rate reduces to the 
Bremermann-Bekenstein bound on information flow [16) . 
Quantum nonequilibrium entropy production. Let us 
consider a closed quantum system, initially in a thermal 
state, whose Hamiltonian Ht is driven by an external 
time-dependent parameter during time r. When this pa- 
rameter is changed in a slow, quasistatic way, the sys- 
tem remains in equilibrium and the process is reversible. 
By contrast, for a fast change, for instance a param- 
eter quench, the system is driven in a nonequilibrium 
state. The nonequilibrium entropy variation associated 
with such a transformation may be defined as [17j . 



ASi, = l3{Wi,) , 



(2) 



where {Wii) — (W) — AF is the difference between the 
total work (W) done on the system during time t and the 
free energy difference AF (the reversible work); as usual 
/3 — l/{kT) denotes the inverse temperature. The irre- 
versible work (VFir) vanishes for a reversible process and 
is defined even if the final state of the system is arbitrar- 
ily far from equilibrium. The probability distribution of 
quantum work is given by the difference EJ^ — E'^ of final 
and initial system energy eigenvalues, averaged over all 
initial states (thermal distribution p° = exp(— /3£'")/Zo) 



and final states (transition probabilities p^ „) [T2], 



An experimental scheme to measure the distribution 
■p (W) in a modulated cold ion trap has been proposed 
in Ref. [IH]- According to Eq. (Is]), the mean work is 
simply (W^) ^ E™,„ (^™ - E°n)Pln..Pl By introduc- 
ing the equilibrium density operator at the final time r, 
p'^^ = exp {—f5Hr)/Zr with eigenvalues p^, we can write 

{W) = 1//3 Y. Pn ^^Pn - 1//3 E P"^™,n ^^^P™ 



l//3 1n(Z,/Zo) 



(4) 



The last term on the right-hand side is equal to AF, 
while the first two are (1//3) times the quantum Kullback- 
Leibler divergence S{pr\\Pr'^), or quantum relative en- 
tropy [in], between the actual density operator of the 
system pr at time r and the corresponding equilibrium 
density operator p'^'^. Using Eq. (pi), we therefore obtain. 



A5i,. = S {pr M"") - tr {pr In pr - Pr In p^'i} 



(5) 



This is an exact expression for the nonequilibrium en- 
tropy production and a quantum generalization of recent 
results presented in Refs. [201 HI]- We note, however, that 
the relative entropy is not a true metric, as it is not sym- 
metric and does not satisfy the triangle inequality; it can 
therefore not be used as a proper quantum distance [22] . 
Furthermore, Eq. ^ is in general difficult to determine 
explicitly. We next derive a lower bound for the quan- 
tum entropy production which we express in terms of the 
fidelity, one of the most commonly used and well-studied 
measures in quantum information theory \lb[ . 

Generalized quantum Clausius inequality. Inequali- 
ties are essential tools of classical and quantum informa- 
tion theory [2^; they allow to express 'impossibilities', 
things that cannot happen, and relate hitherto uncon- 
nected quantities. An elementary example is Klein's in- 
equality, S'(pi||p2) > 0, which asserts the non-negativity 
of the quantum relative entropy jl5) . Combined with 
Eq. ([5]), it immediately leads to the usual Clausius in- 
equality. We shall establish a generalized Clausius in- 
equality by proving that the irreversible entropy vari- 
ation is always larger than the Bures length [T^l E]. 
The Bures metric formally quantifies the infinitesimal 
distance between two density operators as £'^{p+5p, p) = 
iv{5pG}/2, where G obeys pG + G p = 6 p. Dis- 
tances should be physically motivated and to some degree 
unique [J. Wootters' statistical distance, being equal to 
the angle in Hilbcrt space, is the only Riemannian metric 
(up to a constant factor) which is invariant under all uni- 
tary transformations [Hj; it is hence a natural metric on 
the space of pure states. The Bures metric, on the other 
hand, is the generalization of Wootters' metric to mixed 



states [23]; in this sense it represents a natural, unitar- 
ily invariant Riemannian metric on the space of impure 
density matrices [21] . For any two density operators the 
finite Bures length C is given by 



-C (pi, P2) = arccos ( ^jF{pi,p2) 



(6) 



where the fidelity F is defined for an arbitrary pair of 
mixed quantum states as [23 [21] , 



F{pi,P2) 



tr 



Pi P2 \JP\ 



-\ 2 



(7) 



The fidelity is a symmetric, non-negative and unitar- 
ily invariant function, which is equal to one only when 
the two states are identical. For pure quantum states. 
Pi — |'0i)(V'i|i the fidelity reduces to their overlap, 
F{Pi,P2) = tr{piP2} = I (V'ilV'2) P- It has recently been 
shown that if d{pi, P2) is a unitarily invariant norm, then 
the quantum relative entropy satisfies (Ref. [2Z|, Th. 4), 



S{pi\\p2)>'2 



d'^iPi,P2) 



d^t 



1,1 p2,2 



)' 



(8) 



where e*'-' is the matrix with z, j element equal to 1 and 
all other elements 0. Noting that £(e^'^,e^'^) = 7r/2, 
since the two matrices are orthogonal {F{e-^'^, e^'^) = 0), 
we obtain the generalized Clausius inequality, 



AS,,>^C'ipr,p'?) 



(9) 



The quantum entropy production AS'ir is hence bounded 
from below by the geometric distance between the ac- 
tual density operator pr at the end of the process and 
the corresponding equilibrium operator p'^'^, as measured 
by the Bures length; the latter defines a quantum gen- 
eralization of the concept of thermodynamic length [SHE] 
and provides a natural scale to compare AS'ir with. In 
other words, inequality ([9]) quantifies in a precise way 
the intuitive notion that the irreversible entropy produc- 
tion is larger when a system is driven farther away from 
equilibrium. Expression (p| shows that AS^^ is bounded 
by many distances, however, only the Bures length has a 
simple physical interpretation. Equation Q is valid for 
arbitrary quantum processes, including far from equilib- 
rium final states. For infinitesimally close diagonal states, 
we have S{p\\p + dp) ~ 2C'^{p\\p + dp) ~ df{p\\p + dp)/2. 
In the limit of classical, quasi-equilibrium transforma- 
tions, Eq. ^ thus reduces to dS„ > 2/t:'^ dP. 

A simple illustration of the generalized Clausius in- 
equality ([9]) is provided by a time-dependent quantum 
harmonic oscillator, initially at thermal equilibrium. 



Ht = 



2m 



1 2 ■ 
+ -mujf X 



(10) 



and whose angular frequency is varied from wq to wi ac- 
cording to uj'^ = ujI + {ujf - Lol)t/T (cf. Fig. [l| [2H]. 




OJ, 
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FIG. 1: (Color online) Illustration of the generalized Clausius 
inequality (|9| for the quantum harmonic oscillator ( 10 1 with 
linearly varying frequency. Entropy production ASir (red, 
dashed) and lower bound S/tt'^ C^ [pr, p"^) (blue) as a function 
of the final frequency ui\ for h = P = t = 1 and cjq = 0.1. 



This quantum system is exactly solvable and both the 
entropy production and the Bures length can be evalu- 
ated explicitly |29j . The distance from equilibrium is here 
controlled by the ratio wi/wq between final and initial fre- 
quencies. We observe that the relation dSir ~ M"^ /2 is 
valid close to equilibrium, oji/wo — 1, whereas the gen- 
eralized Clausius inequality ([9| holds for all values of wi . 
It is worth noting that for very large Wi , initial and final 
state become maximally distinguishable, i.e. orthogonal, 
and the Bures length (l6| approaches 7r/2. For the pa- 
rameters of Fig. [ll this very far from equilibrium regime 
is reached when the energy of the oscillator is increased 
by more than a factor of twenty. A quantitative estimate 
of the energy change associated with a given irreversible 
entropy production will be given below. 

Quantum entropy production rate. Nonequilibrium 
irreversible phenomena are not only characterized by the 
irreversible entropy change, but also by the rate of en- 
tropy production a — /S.S-„It. The entropy rate a is 
a central quantity that is associated with the speed of 
evolution of a nonequilibrium process [T]. In quantum 
mechanics, the energy of a system imposes a funda- 
mental constraint on its unitary time evolution, as cap- 
tured for instance by the time-energy uncertainty relation 
t > h/AE: the minimal time it takes for a quantum sys- 
tem to evolve to an orthogonal state is always larger than 
the inverse of its initial energy spread AE. A more ac- 
curate expression for the 'quantum speed limit' has been 
derived in Ref. 130! for time-independent Hamiltonians. 
Consider a system initially in state po with mean energy 
Eo - (Ho) and energy spread A^o - i{Hi) - {H^ff/^ 
The minimum time required for the evolution to an ar- 
bitrary state pr is then given by 



max 



2hC^{pr,pa) fi^{pT,Po) 



ttEq 



AEo 



(11) 



where C{pr,po) is the Bures length ^ between pr and 
Pq. The quantum speed limit time is thus entirely de- 
termined by the initial energy (mean or variance) and 
the geometric distance between initial and final states. 
We will now show that Eq. ( [TT] ) sets an upper bound 
to the entropy production rate. We begin by rewriting 
Eq. ^ in the form AS*;, = /3| {Hr) - (Hq) - Fr + Fq\, 
since AS'ir > 0, and focus on the limit of large excitations 
(Ht-) ^ (Hq), which is achieved for long enough driving 
times. By applying the triangle inequality and noting 
that AF = Fr - Fo < (Hr) - (Hq), we obtain 



A5ir < (5{{Hr) + (Ho) + AF) 

< 2f3{{Hr) + (Ho)) ^ 2(3 {Hr) 



(12) 



By combining Eqs. (Illj) and (12), we find that the max- 
imum entropy production rate (Tmax = AS'ir,max/'''i 



IS 



o-max ^ 2/3 (Hr) min 



ttEq AEo 

^2hC'^ iPr,P())' tiC{pr, 



Po) 



(13) 



The above equation expresses the inherent quantum- 
mechanical limit to the nonequilibrium entropy produc- 
tion rate. Equation ( 13 ) simplifies in the limit of far 



from equilibrium transformations when initial and final 
states become orthogonal, £(pr,Po) — 7r/2, and in the 
limit of high temperatures, Eq ~ AEq ~ 1//3, to the 
Bremermann-Bekenstein bound (16j . 



rnr 



(14) 



The Bremermann-Bekenstein bound on information flow 
gives the maximum quantum communication rate (ca- 
pacity) that is possible through a noiseless single channel 
with signals of finite duration. We stress that the present 
derivation is solely based on the thermodynamic defini- 
tion of the entropy production ([2]) and does not make any 
reference to information entropy or channels.; it is thus 
free of the caveats of the original derivations, such as the 
use of the classical Shannon formula for the capacity, or 
the periodic boundary condition approximation |31j . 
The quantum speed limit time (11 1, and hence the en- 



tropy production rate (13), only hold for slowly driven 
systems (i.e. quasi-equilibrium processes), as Eq. (11) 
assumes a time-independent Hamiltonian. For arbi- 



trary nonequilibrium transformations, Eq. (11) can be 
extended using the geometric approach of Ref. |32 



By 

evaluating the time derivative of the angle between initial 
and final state, C {pr, Po), we find that the exact quantum 
speed limit for time-dependent Hamiltonians is 



hC{pr,Po) fl-^iPT;Po) 



Er 



AEr 



(15) 



Contrary to Eq. (11), the minimum time is here de- 



termined by the time averaged mean energy and en- 



ergy 



variance, Er = {1/t) J^ dt {Ht) and AEr 



T 
7.2 r 

4tO 



0.6 




FIG. 2: 

Eq. pi 



(Color online) Quantum speed limit time Tmin, 
(red) and actual driving time r (blue, dashed) for 



the linearly parameterized quantum harmonic oscillator ( 10 1 
with fi = T = 1, 1//3 = and wq = 0.1. 



(l/r) J^ dt {(Hf) - {Htff/^, and not by their initial val- 
ues [33". As a consequence, the quantum speed limit time 
for driven systems can be smaller than for undriven sys- 
tems when Er > Eq {AEr > AEq). This is for instance 



the case at zero temperature: According to Eq. (Ill, a 



quantum system never leaves an initial pure state (infi- 
nite r,„in) in the absence of driving, while Eq. ( 13 ) pre- 
dicts a finite T^in for a driven Hamiltonian. Figure [2] 



shows that for the time-dependent oscillator ( 10 1 at zero 



temperature, the actual driving time r can approach the 
absolute minimum evolution time Tmin within a factor of 
two, for a simple linear change of its angular frequency. 
The general expression for the maximum entropy pro- 
duction rate a for nonequilibrium quantum processes can 



eventually be obtained by combining Eqs. ( 11 ) and (151; 



it is to be regarded as the extension of the Bremermann- 



Bekenstein bound ( 14 ) to arbitrary distances between ini- 



tial and final states, arbitrary initial temperature and 
arbitrary transformation speed. Written in the form 



(Ht) 



> — mm 



Er 



AEr 



hC{pr,Po)' fl'0.{pr,Po) 



(16) 



it provides an estimate for the minimum energy change 
occurring with a given entropy variation in a time r. 

Conclusion. We have developed a generic geometric 
characterization of far from equilibrium quantum pro- 
cesses based on the distinguishability metric on the space 
of quantum states. We have first obtained a generalized, 
more precise form of the Clausius inequality by deriving a 
lower bound for the irreversible entropy production given 
by the Bures length between the nonequilibrium density 
operator of the system and the corresponding equilibrium 
operator. We have further estabhshed the existence of an 
upper bound for the entropy production rate by employ- 
ing the notion of quantum speed limit time which is itself 
a function of the Bures length. The latter is an exten- 



sion of the Bremermann-Bekenstein bound to arbitrary 
nonequilibrium quantum processes. 
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Note that the exact Tmin ( 15 \ does not reduce exactly to 



the estimated Tmin (111 in the limit of slow driving. 



